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Abstract 

Background and Purpose: Treatment outcomes vary greatly in patients with nasopharyngeal carcinoma (NPC). 
The purpose of this study is to evaluate the influence of radiation and chemotherapy drug action pathway gene 
polymorphisms on the survival of patients with locoregionally advanced NPC treated with cisplatin- and fluorouracil- 
based chemoradiotherapy. 

Material and Methods: Four hundred twenty-one consecutive patients with locoregionally advanced NPC were 
prospectively recruited. We utilized a pathway approach and examined 18 polymorphisms in 13 major genes. 
Polymorphisms were detected using the LDR-PCR technique. Multifactor dimensionality reduction (MDR) analysis 
was performed to detect potential gene-gene interaction. 

Results: After adjustment for clinicopathological characteristics, overall survival was significantly decreased in 
patients with the MPO rs2243828 CT/CC genotype (HR=2.453, 95% CI, 1.687-3.566, P<0.001). The ERCC1 
rs3212986 CC (HR=1.711, 95% CI, 1.135-2.579, P=0.010), MDM2 rs2279744 GT/GG (HR=1.743, 95% CI, 
1.086-2.798, P=0.021), MPO rs2243828 CT/CC (HR=3.184, 95% CI, 2.261-4.483, P<0.001) and ABCB1 rs2032582 
AT/AA (HR=1.997, 95% CI, 1.086-3.670, P=0.026) genotypes were associated with poor progression-free survival. 
Prognostic score models based on independent prognostic factors successfully classified patients into low-, 
intermediate-, and high-risk groups. Furthermore, MDR analysis showed no significant interaction between 
polymorphisms. 

Conclusions: Four single nucleotide polymorphisms were associated with survival in patients with locoregionally 
advanced NPC treated with cisplatin- and fluorouracil-based chemoradiotherapy. Combining clinical prognostic 
factors with genetic information was valuable in identifying patients with different risk. 
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Introduction 

Prevalent in the south of China, nasopharyngeal carcinoma 
(NPC) has an ethnic and geographic distribution pattern that is 
distinctive from other head and neck cancers [1,2]. NPC is 
highly responsive to radiotherapy and systemic chemotherapy. 
Based on the results from several large prospective clinical 
trials using cisplatin- and fluorouracil (5-FU)-based 
chemoradiotherapy (CRT) [3-5], the combined treatment has 
become the standard regimen for locoregionally advanced 
NPC. However, even with this multimodal approach, 
approximately 30% of patients experience treatment failure 
within 5 years [4], and treatment outcomes vary widely, even 
among patients of the same clinical stage. Therefore, 
prospectively identifying patients who will have better or worse 
outcomes after CRT would aid in designing appropriate 
treatment strategies, highlighting the need for improved 
predictive markers. 

The expression and activity levels of the critical enzymes 
related to radiation and chemotherapy drug response could 
affect treatment outcomes. Nucleotide excision repair (NER) 
and base excision repair (BER) are two important mechanisms 
involved in the repair of non-specific DNA damage induced by 
radiation and chemotherapy. The major enzymes involved are 
ERCC1 , ERCC2 (NER), and XRCC1 , hOGG1 , APEX1 , ADPRT 
(BER). Common genetic variations of above genes have been 
widely studied in the prognosis and susceptibility of various 
cancers [6-8]. For example, polymorphisms of the XRCC1 
gene, such as rs25487 (Arg399Gln), are capable of altering the 
phenotype of the XRCC1 protein, thus causing a deficiency in 
DNA repair [9] that is significantly related to patient survival [6]. 

Drug metabolism and transportation are important in the 
response of chemotherapy. MTHFR is a key enzyme in the 
folate metabolism pathway that regulates the intracellular folate 
pool for the synthesis and methylation of DNA [10]. High level 
of pre-treatment MTHFR expression was correlated with 
favorable response to fluorouracil-based chemotherapy [11]. 
GSTP1 is a primary enzyme responsible for the detoxification 
of platinum agents. It has been demonstrated to be a predictive 
marker of overall survival (OS) in cancer patients treated with 
cisplatin-based chemotherapy [12]. The ABCB1 gene which 
was also known as MDR1 gene encodes the multi-drug efflux 
pump P-glycoprotein (P-gp) which is involved in the transport of 
a wide range of anti-cancer drugs including cisplatin and 
fluorouracil [13]. The increase expression of ABCB1 was 
related with multidrug resistance and a poor response to 
chemotherapy for increased drugs efflux pathways [14]. 
Polymorphisms of these genes such as MTHFR rs1 801 1 31 , 
GSTP1 rs1 138272 and ABCB1 rs1045642 were associated 
with treatment outcomes in patients receiving chemotherapy 
using 5-Fu and cisplatin [6]. 

VEGF, as a critical angiogenic factor, plays an important role 
in cell growth and survival of endothelial cells and tumor cells. 
Moreover, VEGF production is a potential predictive marker for 
chemotherapy. It was significantly higher in cisplatin-resistant 
cancer cells than that in cisplatin-sensitive parental cells [15]. 
Genetic variation of VEGF was reported to be related with 
clinical outcome of patients treated with platinum-based 



chemotherapy [16]. The FGF-FGFR family also plays a critical 
role in cancer development because of its action in 
angiogenesis [17]. It has been demonstrated that FGFR4 
polymorphism is related with resistance to adjuvant therapy in 
primary breast cancer [18]. 

P53 plays an important role in cellular processes including 
cell-cycle arrest, DNA repair, and apoptotic cell death in 
response to cellular stress including chemotherapy [19]. MDM2 
is a major negative regulator of p53. It directly binds to and 
inhibits p53 by regulating its location, stability, and activity as a 
transcriptional activator [20]. Previous study showed that 
MDM2 polymorphism (rs2279744) located in the first intron of 
the MDM2 promoter was an independent prognostic factor for 
cancer patients [21]. 

MPO is released by neutrophils and macrophages. It is a 
major enzyme involved in generating highly cytotoxic 
hypochlorous acid and other reactive oxygen species (ROS) 
[22], which result in oxidative stress-mediated apoptosis [23]. 
MPO is also extremely important in drug metabolism. With the 
ability to oxidize a wide variety of compounds and a broad 
range of functional groups [24], MPO may enhance the effect 
of chemotherapy. Functional polymorphism which influences 
the expression of MPO could impact on survival in patients 
receiving chemotherapy [25]. 

Because CRTs exert their effects through multistep, 
multigenic cascades, it is unlikely that any single SNP would 
have such a dramatic effect that it could serve as a sole 
predictive marker for response [6]. Therefore, in the current 
study, we use a pathway approach to investigate the impact of 
genetic variations in radiation and chemotherapy drug action 
pathway genes on the survival of patients with locoregionally 
advanced NPC treated with standard CRT. A comprehensive 
panel of 18 SNPs in 13 major genes involved in NER, BER, 
drug metabolism, drug disposition, oxidative stress reaction, 
the p53 pathway, FGFR4 and VEGF were selected. We 
analyzed the association of each polymorphism with patient OS 
and progression-free survival (PFS) individually; then, 
according to multivariate analysis results, we built prognostic 
score models (PSMs) with genotypes and clinical 
characteristics to define different prognostic risk groups. 
Furthermore, the interaction between polymorphisms was also 
investigated. 

Materials and Methods 

Ethic statement 

This study was approved by the independent Institute 
Research Ethics Committee at the Sun Yat-sen University 
Cancer Center (SYSUCC, Guangzhou, P. R. China), and 
written consents were obtained from all participants. 

Patient cohort 

January 2002 and December 2005, 421 consecutive patients 
with locoregionally advanced NPC at the SYSUCC were 
enrolled in the study. All patients received standard CRT 
(details are described below). Patients were considered eligible 
if they had stage III or IVa-b [6th International Union against 
Cancer (UICC)/American Joint Committee on Cancer (AJCC) 
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staging system] histological confirmed NPC. Other eligibility 
criteria included an age of 18 to 65 years, Han Chinese 
ethnicity, an Eastern Cooperative Oncology Group (ECOG) 
performance status of 0 or 1, WHO type II or III NPC, and 
treatment with concurrent CRT using cisplatin with or without 
sequential chemotherapy (induction or adjuvant). Patients who 
had other concomitant malignant diseases or who were 
previously treated with radiotherapy or chemotherapy were 
excluded. 

Treatment and follow-up 

All patients were prospectively included in a disease-specific 
database. The pretreatment evaluation included a complete 
physical examination, magnetic resonance imaging (MRI) of 
the nasopharynx and neck region, fiber optic 
nasopharyngoscopy, chest X-ray, abdominal ultrasound, bone 
scan by emission computed tomography, complete blood count 
and liver and renal biochemistries. 

All patients received definitive radiotherapy with 6-MV 
photons. Uniform radiotherapy protocols for conventional two- 
dimensional radiotherapy (2D-CRT) and intensity-modulated 
radiotherapy (IMRT) at the SYSUCC were followed as 
previously described [26,27]. Radiotherapy was administered 5 
times per week at 2 Gy/d for 2D-CRT and at 2.27 Gy/d for 
IMRT. The accumulated radiation doses were 68-72 Gy to the 
primary tumor, 60-62 Gy to the involved areas of the neck and 
50 Gy to the uninvolved areas. 

Concurrent chemotherapy was administered to all patients. 
Altogether, 28.7% and 25.9% of the patients received induction 
and adjuvant chemotherapy, respectively (Table 1). For 
induction chemotherapy, 2 cycles of PF chemotherapy were 
administered [cisplatin 100 mg/m 2 i.v. drip on day 1 and 5-FU 
1,000 mg/(m 2 d) continuous i.v. for 120 hours, repeated every 3 
weeks]. For concurrent chemotherapy, cisplatin was 
administered at 100 mg/m 2 on days 1, 22, and 43 during 
radiotherapy or at 40 mg/m 2 weekly during radiotherapy. For 
adjuvant chemotherapy, concurrent CRT followed by a 
combination of cisplatin (80 mg/m 2 ) plus 5-FU [1000 mg/(m 2 d)] 
was administered by 96-hour infusion every four weeks for 
three cycles. 

Patients were assessed at the completion of treatment, at 
least once every three months over the first three years and at 
least once every six months thereafter. The patient evaluation 
included a clinical examination, nasopharyngeal endoscopy, 
MRI of the nasopharynx and neck area, chest X-ray and 
abdominal ultrasound. Follow-up ended on September 12, 
2012, with a median follow-up time of 62.0 (range from 5 to 
125) months (Table 1). During the follow-up period, 116 
(27.6%) and 136 (32.3%) patients died or experienced disease 
progression, respectively. Five-year OS and PFS rates for the 
entire patient cohort were 74.3% and 67.8%, respectively. 

DNA extraction and genotyping 

Blood for genotyping was prospectively collected at the time 
of enrollment. All sample collection and storage procedures 
were standardized. Genomic DNA was extracted from 
lymphocytes using the QIAamp DNA Blood Midi Kit (Qiagen, 
Hilden, Germany) according to the manufacturer's protocol. 



Table 1. Demographic and clinical 


characteristics of 


patients. 






Patient characteristics 


N=421 (%) 


Age, years 


Median 


45 


Range 18-65 


<45y 


215 (51.1) 


>45y 


206 (48.9) 


Gender 


Male 


322 (76.5) 


Female 


99 (23.5) 


Histology 


WHO type II 


20 (4.8) 


WHO type III 


401 (95.2) 


T stage a 


T1 


8(1.9) 


T2 


139 (33.0) 


T3 


184 (43.7) 


T4 


90 (21.4) 


N stage a 


NO 


77(18.3) 


N1 


156 (37.1) 


N2 


137 (32.5) 


N3 


51 (12.1) 


Stage a 


III 284 (67.5) 


IVa-b 


137 (32.5) 


EBV DNA 


<4000 copies/ml 


235 (55.8) 


£4000 copies/ml 


186 (44.2) 


Treatment 


CCRT 


191 (45.4) 


IC+CCRT 


121 (28.7) 


CCRT+AC 


109 (25.9) 


RT technique 


2D-CRT 


323 (76.7) 


3D-CRT 


14(3.3) 


IMRT 


84 (20.0) 


Follow-up time (months) 


Median 


62.0 


Range 


5-125 



a 2002 American Joint Committee on Cancer/International Union Against Cancer 
staging system. 



Abbreviation: CCRT, Concurrent Chemoradiotherapy; IC, Induction Chemotherapy; 
AC, Adjuvant Chemotherapy; RT, Radiotherapy; 2D-CRT, Conventional Two- 
Dimensional Radiotherapy; 3D-CRT, Three-Dimensional Conformal Radiotherapy; 
IMRT, Intensity-Modulated Radiotherapy, 
doi: 10.1371/journal.pone.0082750.t001 

The genes involved in drug action pathways were identified 
using the Pharmacogenomics Knowledge Base ( http:// 
www.pharmqkb.org/ ). The genes involved in radiation response 
pathways were identified by a search of the published 
literature. SNPs with a minor allele frequency of 0.05 or more in 
Asian populations were included for evaluation. Overall, we 
selected 18 SNPs from 13 major genes (Table 2). Genotyping 
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of the SNPs was performed by the Shanghai BioWing Applied 
Biotechnology Company (BioWing, Shanghai, China) using the 
ligase detection reaction- polymerase chain reaction (LDR- 
PCR) technique [28,29]. To validate the genotyping results, 
10% of the samples were randomly selected for genotyping by 
a second investigator; the agreement rate was 100%. 

Statistical analysis 

The primary endpoint for this study was OS, defined as the 
time from the date of enrollment to the date of the last follow-up 
visit or death from any cause. The second endpoint was PFS, 
calculated from the date of enrollment to the date of the first 
failure at any site, death from any cause or last follow-up visit. 

Clinicopathological characteristics were dichotomized as 
follows: age (<45 y vs. >45 y), gender (male vs. female), 
histology (WHO type II vs. Ill), T stage (T1-2 vs. T3-4), N stage 
(N0-1 vs. N2-3), clinical stage (III vs. IVa-b), and plasma 
Epstein-Barr virus (EBV) DNA level (<4000 copies/ml vs. 
>4000 copies/ml). Associations of the genotypes with 
clinicopathological characteristics were evaluated by x 2 or 
Fisher's exact tests. The Hardy-Weinberg equilibrium was 
determined for each SNP using a goodness-of-fit x 2 test. The 
impact of the polymorphisms on OS and PFS were examined 
using the Cox proportional hazard model with the calculation of 
hazard ratios (HRs) and 95% confidence intervals (CIs). To 
account for multiple comparisons in the SNP-based analysis, Q 
values set at 0.20 were computed using model-free test (two 
degrees of freedom) P values to quantify the probability that a 
P value may be a false positive [30], accepting a false 
discovery rate (FDR) of 20%. The survival end points were 
analyzed and estimated using the Kaplan-Meier method. The 
significance of the differences among survival curves was 
compared using the log-rank test. Multivariate analyses using 
the Cox proportional hazard model were used to detect 
independent prognostic factors, including genotypes and 
clinicopathological characteristics. The regression coefficient 
("n" in the Cox regression equation HR=e") of each 
independent prognostic factor was then transformed into an 
integral number to build a PSM [31]. We evaluated the 
predictive value of the PSMs and clinical stage by receiver 
operating characteristic (ROC) curve analysis [32]. The areas 
under curves (AUC) were compared between PSMs and 
overall stage. 

To detect potential interactions among polymorphisms, the 
multifactor dimensionality reduction (MDR) analysis was 
performed (MDR software v3.0.2; available on http:// 
sourceforqe.net/proiects/mdr/ ). The non-parametric MDR 
method is described in detail elsewhere [33,34]. Patients with 
missing data for polymorphisms were excluded from the 
analysis. We assumed that patients with beneficial or 
unfavorable genetic profiles have a survival much longer or 
shorter than the median. Therefore, we chose patients in the 
top and bottom quartiles of survival (OS and PFS, respectively) 
in the gene-gene interaction analysis to increase discriminating 
power [35,36]. In the interaction analysis, the ratio between 
patients in the top and bottom survival quartile for each 
genotype combination was evaluated. Combinations with more 
patients in top quartile than in the bottom quartile were 



considered high chance of favorable survival and vice versa. 
This procedure was carried out across 10-fold cross-validation 
samples to avoid over-fitting and was repeated for all possible 
combinations of two to four polymorphisms. The best 
combination was considered, if it had minimal prediction error 
and maximal cross-validation consistency (CVC). Statistical 
significance was further evaluated by a 1000-time permutation 
test to compare observed testing accuracies with those 
expected under the null hypothesis of no association (MDR 
permutation test module v1.0 beta 2; available on http:// 
sourceforqe.net/proiects/mdr/files/mdrpt/ ). 

All statistical tests were two-sided, and a P value of less than 
0.05 was required for statistical significance. All the statistical 
analyses except those specifically mentioned were performed 
using SPSS 16.0 software (SPSS Inc., Chicago, IL). 

Results 

Patient characteristics and distribution of genotypes 

Patient characteristics are presented in Table 1. For each 
SNP, two to seven samples could not be genotyped 
(98.1-99.5% call rate). Because only the C allele was detected 
for rs1 138272 in our patient cohort, it was excluded from the 
analysis. Except for rs3212986 (PO.001), rs1 799793 
(P=0.009), rs2032582 (P=0.020), rs2279744 (P=0.021) and 
rs833061 (P=0.021), all SNPs were tested at the Hardy- 
Weinberg equilibrium (P>0.05). 

Polymorphisms and clinicopathological characteristics 

Table S1 displays the associations between genotypes and 
clinicopathological characteristics. The distributions of XRCC1 
rs25478 and MTHFR rs1801133 were not in equilibration at 
overall stage (P =0.018) and EBV DNA level (P =0.043), 
respectively. For rs25478, genotype GG was significantly more 
frequently observed in patients with stage III (59.2%) vs. stage 
IV disease (47.7%). The rs1801133 genotype CC was more 
frequently observed in patients with high level of plasma EBV 
DNA than in patients with low EBV DNA (51.9% vs. 39.5%, 
respectively). 

Polymorphisms and survival 

Table 2 illustrates the associations between polymorphisms 
and survival in locoregionally advanced NPC patients. 

MPO rs2243828 and MDM2 rs2279744 SNPs were 
significantly associated with patient OS (model-free test, 
PO.001, Figure 1A and P=0.036, Figure 1B, respectively, 
Table 2). Patients with a variant C allele (CT/CC) at rs2243828 
and a variant G allele at rs2279744 (GT/GG) had significantly 
increased death risks relative to those patients with the wild- 
type genotype (rs2243828: HR=2.561, 95% CI, 1.776-3.713, 
P<0.001, Figure 1C; rs2279744: HR=1.778, 95% CI, 
1.074-2.943, P=0.025, Figure 1D). 

We found that the ERCC1 rs3212986, MDM2 rs2279744, 
MPO rs2243828 and ABCB1 rs2032582 SNPs were 
significantly associated with patient PFS (model-free test, 
P=0.023 and Figure 2A, P=0.014 and Figure 2C, PO.001 and 
Figure 2E, P=0.046 and Figure 2G, respectively). Carriers of 
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Table 2. Radiaton and chemotherapy pathway gene polymorphisms and survival in patients with NPC treated with 
cisplatin/5-FU based CRT. 





Pathway 3 


Gene 


SNP 


Genotype 


Overall ; 


survival 






Progression-free survival 














Dead/Alive HR (95%CI) b 


pb 


Q 


Failure/No HR (95%CI) * 


pb 


Q 


R+D 


FGFR4 


rs351855 


























CC 


32/81 


Ref 






40/73 


Ref 












CT 


59/152 


0.952 (0.704-1.288) 


0.751 




72/139 


0.983 (0.668-1 .448) 


0.933 










TT 


21/68 


0.824 (0.558-1.216) 


0.330 




22/67 


0.680 (0.404-1.144) 


0.680 










CT/TT 


80/220 


0.916 (0.686-1.223) 


0.551 




94/206 


0.890 (0.615-1.289) 


0.539 










Model-free 


test c 




0.608 


0.578 






0.273 


0.428 


□+F 


ERCC1 


rs32 12986 


























CC 


78/196 


Ref 






104/170 


Ref 












CA 


27/82 


0.842 (0.542-1 .306) 


0.442 




25/84 


0.567 (0.366-0.878) 


0.011 










AA 


8/24 


1.049 (0.505-2.177) 


0.898 




7/25 


0.597 (0.277-1.284) 


0.186 










CA/AA 


35/106 


0.882 (0.591-1.316) 


0.538 




32/109 


0.573 (0.385-0.853) 


0.006 










Model-free 


test c 




0.721 


0.578 






0.023 


0.105 


D+F 


ERCC2 


rs 1799793 


























GG 


104/263 


Ref 






123/244 


Ref 












AG 


7/37 


0.683 (0.317-1.470) 


0.785 




10/34 


0.727 (0.381-1.387) 


0.333 










AA 


1/4 


1.315 (0.183-9.473) 


0.329 




1/4 


0.823 (0.115-5.901) 


0.846 










AG/AA 


8/41 


0.726 (0.353-1 .494) 


0.384 




11/38 


0.735 (0.396-1.364) 


0.329 










Model-free 


test c 




0.590 


0.578 






0.617 


0.641 


□+F 


ERCC2 


rs13181 


























TT 


93/252 


Ref 






114/231 


Ref 












GT 


17/47 


1.067 (0.636-1.791) 


0.806 




18/46 


0.871 (0.530-1.433) 


0.587 










GG 


1/3 


0.717 (0.100-5.151) 


0.741 




2/2 


1.344 (0.331-5.454) 


0.679 










GT/GG 


18/50 


1.039 (0.627-1.721) 


0.883 




20/48 


0.903 (0.562-1 .453) 


0.675 










Model-free 


test c 




0.915 


0.596 






0.783 


0.715 


R+D+F 


XRCC1 


rs25487 


























GG 


57/173 


Ref 






68/162 


Ref 












AG 


41/110 


1.133 (0.759-1.693) 


0.541 




51/100 


1.158 (0.805-1.665) 


0.428 










AA 


13/20 


1.797 (0.983-3.285) 


0.057 




14/19 


1.587 (0.892-2.823) 


0.116 










AG/AA 


54/130 


1.243 (0.857-1.804) 


0.251 




65/119 


1.230 (0.875-1.728) 


0.234 










Model-free test c 




0.156 


0.576 






0.270 


0.428 


F 


MTHFR 


rs1801 131 


























AA 


55/166 


Ref 






65/156 


Ref 












AC 


49/115 


1.252 (0.851-1.840) 


0.254 




60/104 


1.331 (0.937-1.892) 


0.111 










CC 


8/21 


1.148 (0.547-2.410) 


0.716 




9/20 


1.117 (0.556-2.244) 


0.756 










AC/CC 


57/136 


1.236 (0.853-1.790) 


0.263 




69/124 


1.298 (0.925-1.823) 


0.131 










Model-free 


test c 




0.518 


0.578 






0.280 


0.428 


F 


MTHFR 


rs1801133 


























CC 


47/145 


Ref 






55/137 


Ref 












CT 


52/119 


1.300 (0.876-1.930) 


0.193 




61/110 


1 .293 (0.898-1 .863) 


0.167 










TT 


12/39 


0.948 (0.502-1 .789) 


0.857 




18/33 


1.286 (0.754-2.192) 


0.356 










CT/TT 


64/158 


1.216 (0.834-1.773) 


0.310 




79/143 


1.291 (0.915-1.823) 


0.146 










Model-free 


test c 




0.350 


0.578 






0.348 


0.477 


□+F 


ABCB1 


rs 1045642 


























CC 


44/106 


Ref 






48/102 


Ref 












CT 


51/151 


0.895 (0.597-1 .343) 


0.593 




66/136 


1 .064 (0.734-1 .543) 


0.744 










TT 


19/48 


1.083 (0.631-1.859) 


0.771 




22/45 


1.064 (0.642-1.762) 


0.811 










CT/TT 


70/199 


0.940 (0.643-1.374) 


0.748 




88/181 


1.064 (0.748-1.512) 


0.730 










Model-free test c 




0.742 


0.578 






0.942 


0.782 


D+F 


ABCB1 


rs2032582 


























GG 


30/78 


Ref 






35/73 


Ref 












GT 


45/121 


0.987 (0.621-1.567) 


0.955 




58/108 


1.095 (0.720-1.666) 


0.672 










TT 


20/70 


0.756 (0.430-1.332) 


0.334 




22/68 


0.682 (0.400-1.163) 


0.160 
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Table 2 (continued). 



Pathway 3 


Gene 


SNP 


Genotype 


Overall 


survival 






Progression-free survival 














Dead/Alive HR (95%CI) b 


pb 


Q 


Failure/No HR (95%CI) b 


pb 


Q 








AG 


10/17 


1.279 (0.623-2.623) 


0.502 




8/19 


0.870 (0.403-1.878) 


0.723 










AT 


7/11 


1 .608 (0.704-3.673) 


0.259 




10/8 


2.398 (1.185-4.853) 


0.015 










AA 


0/5 








2/3 


1.150 (0.276-4.785) 


0.848 










GT/TT/AG/AT/AA 


0.953 (0.627-1 .449) 


0.822 






0.996 (0.678-1.464) 


0.984 










Others vs. AT/AA 


1.207 (0.560-2.598) 


0.631 






2.132 (1.176 -3.863) 


0.013 










Model-free 


test c 




0.591 


0.578 






0.046 


0.157 


□ 


MPO 


rs2243828 


























TT 


63/245 


Ref 






71/237 


Ref 












CT 


39/59 


2.131 (1.428-3.178) 


2.09E-4 




53/45 


2.779 (1.946-3.970) 


1.91E-8 










CC 


11/1 


8.994(4.720-17.136) 


2.42E-11 




12/0 


12.041 (6.457-22.454) 


5.03E-15 










CT/CC 


50/60 


2.561 (1.776-3.713) 


7.14E-7 




65/45 


3.226 (2.302-4.521) 


1.01E-11 










Model-free 


test c 




1.88E-11 


2.08E-10 






4.03E-17 


5.52E-16 


R 


hOGG1 


rs 1052 133 


























GG 


39/106 


Ref 






44/101 


Ref 












CG 


57/148 


1.045 (0.695-1.572) 


0.832 




76/129 


1.257 (0.886-1.825) 


0.228 










CC 


16/49 


0.899 (0.502-1 .609) 


0.720 




15/50 


0.725 (0.404-1.304) 


0.283 










CG/CC 


73/197 


1 .009 (0.683-1 .490) 


0.964 




91/179 


1.121 (0.781-1.608) 


0.536 










Model-free 


test c 




0.867 


0.596 






0.114 


0.312 


R 


APEX1 


rs 1130409 


























TT 


47/107 


Ref 






52/102 


Ref 












GT 


51/151 


0.824 (0.554-1 .225) 


0.338 




64/138 


0.962 (0.667-1.387) 


0.835 










GG 


14/44 


0.805 (0.443-1 .462) 


0.476 




19/39 


1.011 (0.597-1.709) 


0.969 










GT/GG 


65/195 


0.820 (0.563-1.193) 


0.278 




83/177 


0.973 (0.688-1.376) 


0.875 










Model-free 


test c 




0.582 


0.578 






0.970 


0.782 


R 


ADPRT 


rs1 136410 


























TT 


40/89 


Ref 






46/83 


Ref 












CT 


51/156 


0.796 (0.526-1 .204) 


0.280 




63/144 


0.834 (0.570-1.220) 


0.349 










CC 


22/56 


0.928 (0.552-1.562) 


0.779 




27/51 


1.023 (0.636-1.646) 


0.926 










CC/CT 


73/212 


0.832 (0.565-1 .223) 


0.349 




90/195 


0.883 (0.618-1.260) 


0.491 










Model-free 


test c 




0.544 


0.578 






0.541 


0.618 


D 


MDM2 


rs2279744 


























TT 


18/80 


Ref 






21/77 


Ref 












GT 


49/134 


1 .586 (0.923-2.722) 


0.095 




61/122 


1.735 (1.056-2.851) 


0.030 










GG 


46/89 


2.041 (1.183-3.521) 


0.010 




54/81 


2.127 (1.284-3.524) 


0.003 










GT/GG 


95/223 


1.778 (1.074-2.943) 


0.025 




115/203 


1.899 (1.192-3.026) 


0.007 










Model-free 


test c 




0.036 


0.199 






0.014 


0.096 


D 


VEGF 


rs20 10963 


























GG 


45/113 


Ref 






48/110 


Ref 












CG 


47/146 


0.877 (0.583-1.321) 


0.530 




61/132 


1.069 (0.732-1.561) 


0.729 










CC 


20/45 


1.084(0.640-1.836) 


0.764 




25/40 


1 .254 (0.773-2.035) 


0.359 










CG/CC 


67/191 


0.930 (0.637-1 .358) 


0.708 




86/172 


1.117 (0.784-1.591) 


0.539 










Model-free 


test c 




0.687 


0.578 






0.655 


0.641 


D 


VEGF 


rs833061 


























TT 


58/167 


Ref 






70/155 


Ref 












CT 


49/124 


1.144 (0.782-1.673) 


0.489 




60/113 


1.159 (0.821-1.637) 


0.401 










CC 


5/12 


1.119 (0.448-2.794) 


0.809 




4/13 


0.721 (0.263-1.975) 


0.524 










CT/CC 


54/136 


1.141 (0.788-1.673) 


0.485 




64/126 


1.117 (0.796-1.568) 


0.523 










Model-free 


test c 




0.783 


0.578 






0.519 


0.618 


D 


VEGF 


rs3025039 


























CC 


79/212 


Ref 






100/191 


Ref 












CT 


32/83 


0.951 (0.630-1.435) 


0.811 




31/84 


0.721 (0.481-1.079) 


0.111 










TT 


1/7 


0.387 (0.054-2.780) 


0.345 




3/5 


0.957 (0.303-3.021) 


0.941 










CT/TT 


33/90 


0.883 (0.586-1.332) 


0.553 




34/89 


0.482 (0.644-1 .060) 


0.095 
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Table 2 (continued). 





Pathway 3 Gene 


SNP 


Genotype Overall survival 






Progression-free survival 










Dead/Alive HR (95%CI) 6 


P b 


Q 


Failure/No HR (95%CI) 6 


P b 


Q 






Model-free test c 


0.630 


0.578 




0.281 


0.428 



a R, Radiation; D, DDP; F, 5-Fluorouracil 

b HR, 95%CI and P values were calculated by using Cox proportional hazard model with no adjustment. 
c The association analysis was done based on the genetic model-free test (2 degrees of freedom) 
doi: 10.1371/journal.pone.0082750.t002 




Figure 1. Kaplan-Meier curves for overall survival in patients with locoregionally advanced nasopharyngeal carcinoma 
treated with cisplatin- and fluorouracil-based chemoradiotherapy. (A) Survival time based on the MPO rs2243828 TT, CT and 
CC genotypes. (B) Survival time based on the MDM2 rs2279744 TT, GT and GG genotypes. (C) Survival time based on the MPO 
rs2243828 TT and CT/CC genotypes. (D) Survival time based on the MDM2 rs2279744 TT and GT/GG genotypes. 

doi: 10.1371/journal.pone.0082750.g001 



the ERCC1 rs3212986 (C8092A) genotypes CA/AA had 
significantly improved PFS compared with the wild-type 
homozygote CC (HR=0.573, 95% CI, 0.385-0.853, P=0.006, 
Figure 2B). Patients with the MDM2 rs2279744 (SNP209TG) 
genotype GT/GG had a 1 .8-fold greater progression risk than 
those patients with the TT genotype (HR=1.899, 95% CI, 
1.192-3.026, P=0.007, Figure 2D). Compared to carriers of the 
MPO rs2243828 (T-764C) TT genotype, the 5-year PFS rate 
was significantly lower for CT/CC carriers (HR=3.226, 95% CI, 
2.302-4.521, P<0.001, Figure 2F). Patients with the ABCB1 
rs2032582 AT/AA genotype had significantly increased 



progression risk relative to patients with other genotypes 
(GG/GT/TT/AG vs. AT/AA, HR=2.132, 95% CI, 1.176-3.863, 
P=0.013, Figure 2H). 

Multivariate analysis 

Multivariate analysis identified that age>45 y (HR=1.898, 
95% CI, 1.293-2.785, P=0.001), N2-3 (HR=1.593, 95% CI, 
1.068-2.376, P=0.022), EBV DNA>4000 copies/ml (HR=1.579, 
95% CI, 1.075-2.320, P=0.020) and MPO rs2243828 CT/CC 
genotype (HR=2.453, 95% CI, 1.687-3.566, P<0.001) were 
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Figure 2. Kaplan-Meier curves for progression-free survival in patients with locoregionally advanced nasopharyngeal 
carcinoma treated with cisplatin- and fluorouracil-based chemoradiotherapy. (A) Survival time based on the ERCC1 
rs3212986 CC, CA and AA genotypes. (B) Survival time based on the ERCC1 rs3212986 CC and CA/AA genotypes. (C) Survival 
time based on the MDM2 rs2279744 TT, GT and GG genotypes. (D) Survival time based on the MDM2 rs2279744 TT and GT/GG 
genotypes. (E) Survival time based on the MPO rs2243828 TT, CT and CC genotypes. (F) Survival time based on the MPO 
rs2243828 TT and CT/CC genotypes. (G) Survival time based on the ABCB1 rs2032582 GG, GT, TT, AG, AT and AA genotypes. 
(H) Survival time based on the ABCB1 rs2032582 GG/GT/TT/AG and AT/AA genotypes. 

doi: 10.1371/joumal.pone.0082750.g002 
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Table 3. Multivariate analysis of treatment outcomes. 



Variables 


HR 


95% CI 


P 


Overall survival 


Age (<45 vs. >45y) 


1.898 


1.293-2.785 


0.001 


Gender (Male vs. Female) 


0.766 


0.479-1 .224 


0.264 


Histology (WHO type II vs. WHO type III) 


1.441 


0.581-3.574 


0.431 


T stage (T1-2VS. T3-4) 


1.238 


0.832-1 .840 


0.292 


N stage (NO-1 vs. N2-3) 


1.593 


1.068-2.376 


0.022 


EBV DNA level (<4000 vs. >4000 copies/ml) 


1.579 


1.075-2.320 


0.020 


MPO rs2243828 (TT vs. CT/CC) 


2.453 


1.687-3.566 


<0.001 


MDM2 rs2279744 (TT vs. GT/GG) 


1.550 


0.930-2.582 


0.092 


Progression-free survival 


Age (<45 vs. >45y) 


1.040 


0.738-1.464 


0.823 


Gender (Male vs. Female) 


0.754 


0.495-1.149 


0.189 


Histology (WHO type II vs. WHO type III) 


0.944 


0.456-1 .955 


0.878 


T stage (T1-2 vs. T3-4) 


1.204 


0.840-1.725 


0.313 


N stage (NO-1 vs. N2-3) 


1.415 


0.980-2.045 


0.064 


EBV DNA level (<4000 vs. >4000 copies/ml) 


1.512 


1.059-2.160 


0.023 


MPO rs2243828 (TT vs. CT/CC) 


3.184 


2.261-4.483 


<0.001 


ERCC1 rs3212986 (CA/AAvs. CC) 


1.711 


1.135-2.579 


0.010 


MDM2 rs2279744 (TT vs. GT/GG) 


1.743 


1.086-2.798 


0.021 


ABCB1 rs2032582 (others vs. AT/AA) 


1.997 


1.086-3.670 


0.026 



doi: 10.1371/journal.pone.0082750.t003 



negative independent factors for OS (Table 3). This analysis 
also indicated that EBV DNA>4000 copies/ml (HR=1.512, 95% 
CI, 1.059-2.160, P=0.023), ERCC1 rs3212986 CC genotype 
(HR=1.711, 95% CI, 1.135-2.579, P=0.010), MDM2 rs2279744 
GT/GG genotype (HR=1.743, 95% CI, 1.086-2.798, P=0.021), 
MPO rs2243828 CT/CC genotype (HR=3.184, 95% CI, 
2.261-4.483, P<0.001) and ABCB1 rs2032582 AT/AA genotype 
(HR=1.997, 95% CI, 1.086-3.670, P=0.026) were negative 
independent factors for PFS (Table 3). 

Prognostic score models 

To build systemic PSMs for OS and PFS, an integral score 
was derived from the regression coefficients of each 
independent prognostic factor. If the factor was absent, a score 
of zero was recorded. If the factor was present, a score of 1 to 
3 was recorded according to the n value (Table 4). The 
maximum scores for OS and PFS were 5 and 8, respectively. 
The overall score of the PSM for each patient was calculated 
as the total of the scores of each independent factor. All 
patients were then categorized into three groups according to 
cut-off points at the 25 th and 75 th percentiles of the score 
distribution as follows: low-risk group (total score 0-1 for OS, 
0-2 for PFS), intermediate-risk group (total score 2-3 for OS, 
3-4 for PFS) and high-risk group (total score 4-5 for OS, 5-8 for 
PFS). For OS, 235, 102 and 84 patients were in the low-, 
intermediate- and high-risk groups, respectively, with 5-year 
OS rates of 88.0%, 71.8% and 37.2%, respectively (PO.001, 
Figure 3A). Referring to PFS, 176, 193 and 52 patients were in 
the low-, intermediate- and high-risk groups, respectively, and 
the 5-year PFS rates were 81.8%, 58.9% and 40.4%, 
respectively (P<0.001, Figure 3B). Figures 4A and 4B show the 
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Table 4. Prognostic score models. 





Overall survival 






Progression-free survival 




= actors 


n (HR=< 


; n )ScoreFactors 


n (HR=< 


5 n )Score 


Age (>45y) 


0.64 


1 


EBV DNA (>4000 
copies/ml) 


0.41 


1 


N stage (N2-3) 


0.47 


1 


MPO rs2243828 
(CT/CC) 


1.16 


3 


EBV DNA (>4000 
copies/ml) 


0.46 


1 


ERCC1 rs3212986 (CC) 


0.54 


1 


MPO rs2243828 
(CT/CC) 


0.90 


2 


MDM2 rs2279744 
(GT/GG) 


0.56 


1 








ABCB1 rs2032582 
(AT/AA) 


0.69 


2 



doi: 10.1371/journal.pone.0082750.t004 



ROC curves for the PSMs and overall stage. The AUCs of the 
PSMs for OS and PFS were significantly larger than the AUCs 
for overall stage (OS: 0.702 vs. 0.572, respectively, P=0.002; 
PFS: 0.745 vs. 0.571, respectively, P<0.001) 

Multifactor dimensionality reduction analysis 

We further performed the MDR analysis to evaluate potential 
interactions between polymorphisms. As shown in Table 5, the 
overall best MDR model was the one factor model using MPO 
rs2243828, which had a maximum CVC (10/10 for both OS and 
PFS) and minimum prediction error (43% for OS and 36% for 
PFS). Moreover, no potential interactions were found between 
investigated polymorphisms in this study (permutation test 
P>0.05 for the best two to four polymorphisms combinations). 

Discussion 

In this study, we used a pathway approach to systematically 
investigate the associations between the genotypes of 18 
SNPs in 13 major genes from radiation and chemotherapy 
response pathways and survival in patients with locoregionally 
advanced NPC treated with standard CRT. Following 
adjustment for clinicopathological characteristics, we found that 
the MPO rs2243828 SNP was independently significant 
associated with OS and PFS and that the ERCC1 rs3212986, 
MDM2 rs2279744, and ABCB1 rs2032582 SNPs were 
independently significant associated with PFS. Moreover, we 
performed MDR analysis and found no significant interaction 
between polymorphisms. We believe our results to be 
biologically plausible for the following evidences. 

As mentioned above, MPO could oxidize a wide variety of 
compounds and a broad range of functional groups [24], then 
enhance the effect of chemotherapy. MPO rs2243828 (T-764C) 
is a functional promoter polymorphism that is linked with 
rs2333227 (G-463A) [6]. The wild-type alleles (T and G) are 
associated with higher expression levels of MPO [37,38]. 
Previous studies have found that low-activity MPO genotypes 
are associated with poor survival in esophageal and breast 
cancer patients treated with chemotherapy [6,25]. We also 
showed that patients with a variant allele had approximately 
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Figure 3. Kaplan-Meier survival curves for the three different risk groups classified by prognostic score models. (A) 

Overall survival. (B) Progression-free survival. 

doi: 10.1371/journal.pone.0082750.g003 



2.5-fold greater death risks and 3.2-fold greater progression 
risks than wild-type homozygous patients. To the best of our 
knowledge, this is the first study focusing on the relationship 
between the polymorphisms of MPO and NPC prognosis. 



However, some studies have found that the variant allele has a 
protective effect against cancer [39,40]. We hypothesize that 
this contradictory finding is a result of the variations in ROS 
activity during the different phases of cancer. In the early 
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1 - Specificity 
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Figure 4. Receiver operating characteristic curves of prognostic score models (PSMs) and overall stage. (A) PSM and 

overall stage for overall survival. (B) PSM and overall stage for progression-free survival. 

doi: 10.1371/joumal.pone.0082750.g004 
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Table 5. MDR analysis for the survival prediction in patients 
with locoregionally advanced nasopharyngeal carcinoma. 





The best combination 


TBA 


cvc 


P value a 


Overall survival 


rs2243828 


0.57 


10/10 


0.30 


rs2243828, rs2032582 


0.56 


8/10 


0.42 


rs2032582, rs2279744, rs1 136410 


0.44 


3/10 


nr" 


rs2032582, rs1136410, rs1801131, rs2010963 


0.43 


3/10 


nr" 


Progression-free survival 


rs2243828 


0.64 


10/10 


0.06 


rs2243828, rs1801 131 


0.62 


6/10 


0.11 


rs2243828, rs1801131, rs2032582 


0.53 


6/10 


0.67 


rs2243828, rs2032582, rs2279744, rs351855 


0.45 


4/10 


nr" 



The above SNPs were specified as following: MPO rs2243828, ABCB1 rs2032582, 
MDM2 rs2279744, ADPRT rs1136410, VEGF rs2010963, MTHFR rs1801131, 
FGFR4 rs351855. 

P value for 1000-fold permutation test. 3 

b. Negative result, P value nearly equals to 1.00. 

Abbreviation: MDR, Multifactor dimensionality reduction; TBA, Testing balanced 
accuracy; CVC, Cross-validation consistency, 
doi: 10.1371/journal.pone.0082750.t005 

phases of carcinogenesis, ROS are important molecules that 
are involved in damaging normal cells and thus contribute to 
the induction of carcinogenesis. On the contrary, in the context 
of cancer treatment, ROS play a different role by contributing to 
tumor cell death mechanisms. 

High expression of MDM2 has been detected in several 
cancers and is related to decreased response to treatment and 
poor prognosis [41]. MDM2 rs2279744 (SNP309 T>G) is 
located in the first intron of the MDM2 promoter and therefore 
has the potential to influence the expression of MDM2. 
Previous studies have shown that cells carrying the GG 
genotype have an increased binding affinity for the 
transcriptional activator Sp1, which subsequently results in 
higher expression levels of MDM2 mRNA and protein and 
elicits a decreased response to DNA-damaging agents [42]. 
Several studies have also demonstrated that the GG genotype 
in patients correlates with poor survival [21]. Liu et al found that 
MDM2 T309G (rs2279744) genotypes were not related with OS 
and PFS in patients with advanced non-small cell lung cancer 
(NSCLC) treated with platinum-based chemotherapy, however, 
the combined analysis showed that significant shorter survival 
was in patients with the p53 Pro/Pro and MDM2 GG genotype 
[43]. With respect to NPC, rs2279744 is associated with the 
risk of developing the disease [44]; however, no studies have 
focused on the clinical outcome of NPC related to this SNP. In 
the present study, we found that patients with the variant allele 
G had significant lower PFS than those patients with the wild- 
type homozygote. This association between genotypes and 
PFS could be explained by the attenuated activity of p53; this 
attenuation is related to the poor chemo/radio-sensitivity of the 
tumor [45,46]. 

The NER pathway is important in the repair of the DNA 
adducts which were typically caused by cytotoxic drugs such 
as cisplatin. ERCC1 is a major component of the NER complex 



and is the rate-limiting enzyme in the pathway. Previous 
studies have found that over expression of ERCC1 is 
associated with resistance to cisplatin-based chemotherapy in 
various cancers [47], including NPC [48]. A common 
polymorphism of the ERCC1 gene, rs3212986 (C8092A) is 
located in the 3' untranslated region and therefore may affect 
mRNA stability and result in a decreased expression level [49]. 
However, the association of the rs32 12986 genotypes and 
cancer prognosis has been inconsistent. In the present study, 
we found that patients with the wild-type homozygous allele 
had significantly lower PFS rates than did those patients with 
CA/AA genotypes. Our result is consistent with study using 
cisplatin-based treatment for gastric cancers [50]. However, 
other studies have reported the opposite results in NSCLC and 
ovarian cancer patients who were treated with platinum-based 
chemotherapy [51,52]. These contradictory results may be 
ascribed to differences in patient populations, tumor 
characteristics or treatment combinations. Altough ERCC2 is 
another major component of NER, we did not find any 
association between ERCC2 genotypes and survival in NPC 
patients. A previous study also found ERCC2 polymorphism 
was not related with survival in NSCLC patients receiving 
platinum-based chemotherapy [53]. 

The ABCB1 is closely related to clinical multidrug resistance. 
The mechanism has multiple aspects, including reduced drug 
accumulation, altered drug metabolism, increased tolerance of 
cellular damage and diminished apoptotic signaling [54]. The 
rs2032582 (G>T/A) SNP is located in exon 21 of the ABCB1 
gene and generates an amino acid change from Ala to Thr/Ser. 
This SNP is often at linkage disequilibrium with another SNP at 
exon 26, rs1 045642 (3435 C/T) [55]. The phenotypes of both 
rs2032582 and rs1 045642 are associated with the clinical 
outcomes of several cancers [56,57]. In the present study, we 
did not find any correlation between the rs1 045642 genotypes 
and patient survival but did find that patients harboring the 
rs2032582 AT/AA genotypes had 2.13-fold greater progression 
risks than other genotypes. Although the exact biological 
reason is unclear, the change in activity of P-gp may account 
for the change in chemotherapy sensitivity. 

Additionally, we found that the XRCC1 rs25478 and MTHFR 
rs1801133 genotypes were correlated with clinical stage and 
EBV DNA level, respectively. As mentioned above, XRCC1 is a 
key protein that is directly involved in the repair of DNA base 
damage, and the Arg399Gln amino acid variant may alter the 
phenotype of the XRCC1 protein, causing deficient DNA repair 
[9]. XRCC1 Arg399Gln was reported to be associated with 
susceptibility and prognosis of various cancers [6, 8]. The 
rs1801133 (C667T) variant is a common polymorphism of the 
MTHFR gene that leads to an amino acid substitution and 
decreased enzyme activity [58]. A previous study showed that 
patients with the variant allele CT/TT have improved survival 
compared with patients with the CC allele [59]. However, the 
reasons for the association of these two SNPs with the 
identified NPC characteristics are still unclear and require 
intensive investigation. 

Liu et al performed a multi-loci analysis in NSCLC patients 
and they found that interactions among XRCC1 Arg194Trp, 
XPC PAT, FAS G-1377A, and FASL T-844C were associated 
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with sensitivity to platinum-based chemotherapy [60]. However, 
the overall best MDR model was the one factor model using 
MPO rs2243828 and permutation test P values for the best two 
to four polymorphisms combinations were over 0.05 in the 
present study, revealing no significant interactions among 
polymorphisms. 

PSMs were built for OS and PFS containing the significant 
SNPs described above. These PSMs categorized patients into 
three risk groups according to their prognosis and 
demonstrated certain value in predicting survival. 

There are some limitations to this study, however. First, as 
with any study of modest size, this study may lack a degree of 
generalizability. Second, 18 SNPs were examined in our study, 
which could lead to false results due to multiple comparisons; 
to address this concern, we used a Q value to maintain the 
FDR under 20%. Finally, our results require validation with a 
large patient cohort prior to clinical application. Therefore, an 
additional prospective multicenter study should be conducted 
to further validate our results within the NPC patient population. 

In conclusion, this study used a pathway approach to 
demonstrate that genetic variations within MPO, MDM2, 
ERCC1 and ABCB1 were associated with survival in patients 
with locoregionally advanced NPC treated with cisplatin- and 5- 
FU-based CRT. Furthermore, our PSMs demonstrated that 
genetic polymorphisms in combination with clinical prognostic 
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